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We consider scattering state contributions to the partition function of a two-dimensional (2D) 
plasma in addition to the bound-state sum. A partition function continuity requirement is used 
to provide a statistical mechanical heuristic proof of Levinson's theorem in two dimensions. We 
show that a proper account of scattering eliminates singularities in thermodynamic properties of the 
nonideal 2D gas caused by the emergence of additional bound states as the strength of an attractive 
^vj^ ' potential is increased. The bound-state contribution to the partition function of the 2D gas, with a 

weak short-range attraction between its particles, is found to vanish logarithmically as the binding 
energy decreases. A consistent treatment of bound and scattering states in a screened Coulomb 
' potential allowed us to calculate the quantum-mechanical second virial coefficient of the dilute 2D 

electron-hole plasma and to establish the difference between the nearly ideal electron-hole gas in 
GaAs and the strongly correlated exciton/free-carrier plasma in wide-gap semiconductors such as 
ZnSe or GaN. 
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I. INTRODUCTION 



Two-dimensional (2D) systems play a central role in contemporary condensed matter physics. Novel phenomena 
such as the quantum Hall effect [|| observed when a 2D electron gas at low temperature is subjected to a strong 
^ magnetic field, as well as practical developments based on quasi-2D systems, e.g., high-mobility field-effect transistors 
1^ or semiconductor quantum-well lasers Q| have brought significant technological advances. Such devices are based 
on the quasiequilibrium response of the internal electron or electron-hole plasmas to an external stimulation. Hence it 
is essential to understand the fundamental quantum-statistical properties of such two-dimensional interacting plasmas 
, at finite temperatures. 

■ One of the well-known differences between 2D and 3D nonrclativistic quantum mechanics is the presence of at least 
CO one bound state for any symmetric attractive potential in two dimensions. This bound state, with binding energy 
Eh, brings a non-vanishing contribution of exp(Efi/kBT) to the two-body part of the partition function even if 
the interaction is weak and the state is very shallow. For a dilute gas this contribution introduces a deviation from 
the ideal gas law that is larger than the correction due to the Fermi or Bose statistics of the particles. However, 
it is clear that 2D gases with a vanishing inter-particle interaction strength should be well described by free Fermi 
or Bose gas models. This contradiction will be resolved in the present paper. A related question is what happens 
to the partition function when additional bound states appear with increasing strength of interaction between the 
, particles. In the 3D case the answer to this question is based on a careful consideration of states in the continuum, 
'"^ ■ which are modified by the interaction, i.e., taking scattering into account in the partition function calculation [Q,^. 
^ ] The same approach should be applied in two dimensions, however, scattering theory in two dimensions is relatively 
O undeveloped compared to the 3D case. For example, the relation between low-energy scattering and bound states, 
which has important consequences in the statistical mechanics, has been considered only very recently in two 
dimensions. In the present paper the connection between scattering and the statistical mechanics of a 2D plasma is 
studied. 

In the next section we introduce the relation of the two-dimensional scattering phase shift to the partition function 



• and show that a proper consideration of the scattering states removes discontinuities in the partition function in an 



analogous fashion to the 3D case. In Sects. HI and IV we explore this in more detail using an analytical model with an 
attractive square well potential before turning to a more realistic model of the statically screened electron-hole plasma, 
which is the main focus of this paper. We also compare in both cases the influence on the second virial coefficient 
of the exchange interaction and the screened direct interaction. Such considerations are crucial in understanding the 
nature of the strongly correlated electron- hole plasma in semiconductor quantum wells. 
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II. PARTITION FUNCTIONS AND LEVINSON'S THEOREM 



The two-body interaction part of the partition function of 2D interacting Bohzmann particles is given by 

Zint = y^exp(-£',n^^/fc_BT) 



1 r g exp(-,V4) dg, (1) 



^ rvvfc rv^ 



m— — oo 



where = 2iikBT/h^ , is the reduced mass, m is the projection of the angular momentum onto the axis 
normal to the plane of 2D motion (m = 0, ±1, ±2, . . .), Smiq) is the 2D scattering phase shift Q dependent on the 
relative-motion momentum q, Em,v are the bound-state energies (index u enumerates bound states with given to), 
and the double sum in the first term ranges only over bound states. Eq. (|^) is the 2D analogue of the Beth-Uhlenbeck 
formula and can be derived in the same fashion as in the 3D case ^ . 

Often only the first term in Eq. (Q) is considered when calculating an internal partition function, neglecting the 
phase shift term. For an attractive potential gU{p) , as g decreases bound state energies increase towards the 
continuum. As such a state reaches zero energy, a partition function that contains only the bound state sum will be 
discontinuous. These unphysical singularities would extend to all the thermal properties, such as pressure and specific 
heat. Integrating by parts we can rewrite Eq. as 



y \^Z^^rn{q) \ cxp(-g2/gf|,)qd(7 . 



+ / I > ^ I cM-q'lqT)qdq . (2) 



For nonzero temperature the integral term in the right-hand side of Eq. (^) is a smooth function of the interaction 
strength g. If the phase shifts satisfy the condition: 

lim 8ra{q) = VmT^ , (3) 

where Vm is the number of bound states with given to , the zero-energy part of the phase shift integral in Eq. (|^) 
exactly cancels the zero-energy part of the bound state sum, removing the discontinuity in Zint as a function of 
the interaction strength. This cancellation is similar to the well known behaviour in 3D where the partition function 
discontinuities are removed |^ with the help of Levinson's theorem Equation (||) constitutes the 2D statement 
of Levinson's theorem. 

As a central theorem of scattering theory [ pT| , Levinson's theorem has been discussed for Dirac particles, multichan- 
nel scattering, multi-particle single-channel scattering, one-dimensional scattering systems, impurities in Aharonov- 
Bohm rings, systems with non-uniform effective mass, and even for time-periodic potentials p^ . However, its appli- 
cability to the 2D scattering problem has been considered only recently. In Ref. |^ the 2D statement of Levinson's 
theorem, Eq. (|^), was proposed and verified empirically, while in Ref. [Q this theorem was more rigorously estab- 
lished for cutoff potentials using the Green-function method. The above arguments, based on the partition function 
continuity requirement, provide an additional statistical mechanical justification of Levinson's theorem in 2D. 

The two-body interaction part of the partition function can be used to calculate the second virial coefficient, B{T), 
which characterizes the first correction to the ideal gas law in the low-density expansion of pressure 

P = nksT {1 + Bn + . . .) . (4) 

B is positive for repulsive potentials, causing an increase in the pressure over its ideal-gas value, and negative for 
attractive potentials, causing a decrease in the pressure. A calculation of the second virial coefficient is meaningful 
in the dilute gas regime, where the mean inter-particle spacing is larger than the thermal wavelength A = 

{2ttTt? /MkBTY^"^ and the higher order terms in Eq. are negligible. For a free, Bose or Fermi gas in two dimensions, 
B{T) = =fA^/4 the plus sign applying to the Fermi case, as the Pauli principle introduces an effective repulsion 
between fermions, thereby increasing the pressure. There is extensive ongoing research in the statistical mechanics 
of anyons, 2D particles obeying fractional statistics ||lj], and the second virial coefficient of a free anyon gas lies in 
between the bosonic and fermionic value [ p^ . In this paper we focus on the relative importance of interaction-induced 
bound and scattering states on the second virial coefficient. 
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The 3D analysis of the second virial coefRcient is easily reformulated for a 2D interacting gas . For a system 
of identical particles with spin s the second virial coefficient is 



A2 



2s + 1 



1 



2Z, 



2s + 1 



(5) 



where the upper sign is for bosons and the lower sign for fermions. The exclusion principle modifies the sum over m 
in Eq. depending on the angular momentum parity, and the partition function (using Eq. (|^)) is 



1±(-1)' 



2 f°° 
+ (2s + 1) — 5- / qdq 

t^Qt Jo 



E 



_rn— — oo 



[eM-Em,./kBT) - 1] 

i±(~ir\ ■ 

Om(q) 



exp(-g /qt) 



(6) 



The electron-hole plasma constitutes a mixture of two components and for a binary mixture of components C 
ad D having si 
coefficient is fl^ 



(s) (s') 

and D having second virial coefRcients Bq and Bj^ and densities nc and nu, respectively, the second virial 



BiT) = (^YbI^^ + 2 



\ n / 



ncriD 



B. 



CD 



D 



In Eq. (0) n = nc + njj and 



BcD — —^^Zint 



where 



2 ^ 2Trh 



McM^ 



D 



Mc + Md ' 

and Zint is given by Eq. (Ill) or Eqs. (M and (|3|) with the properly chosen reduced mass /i. 



(7) 



(8) 
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III. BOLTZMANN GAS WITH SHORT-RANGE ATTRACTION 



Our first example, the Boltzmann gas with weak short-range attraction, is chosen to elucidate how the second virial 
coefRcient at given temperature vanishes when the binding energy Eb decreases, even though the bound-state part of 
the internal partition function Zbound = exTp^Eb/ksT) approaches unity rather than vanishes. To trace the precise 
nature of the cancellation of Zbound it is convenient to use Eq. (|l|) for Zint without the application of Levinson's 
theorem. 

Let us assume that the 2D particles interact via an attractive square- well potential of radius a and depth Vq. 
This simple model allows analytical treatment which provides insight into the generic behaviour of a gas of attracting 
particles. To evaluate the partition function we need to analyse both the bound and scattering states in this potential. 
The binding energies for any value of angular momentum m can be easily found by matching the logarithmic derivative 
of the radial wave function at p — a : 

-k2J|„.| + i (a^K^-K^^ ^ KK\„,\ + l{Ka) ^^^^ 



\m\ 



(a 74^) ^|r„|(Ka) 



where Kq = 2fiVo/Ti , = 2fiEb/?i , Jm{x) is the Bessel function of the Rrst kind and Km{x) is the modified 
Bessel function of the second kind. Note that for < 2.4 (i.e. smaller than the first root of Jo{x) ) there is only 
one bound state (having m = 0) and for < 1 this state is very shallow, e.g. for kqa = 1 , Eb/Vo ~ 0.04 , for 
KqU = 0.5 , Eb/Vo w 2 X 10^^. For a shallow m — state the transcendental equation for binding energies Eq. ( [lO| ) 
reduces to 

KQaJi{KQa) KaKi^Ka) 1 
Jo{Koa) KqIho) In(cKa) ' 
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where c = exp(7)/2 ( 7 w 0.5772157 .. . is Euler's constant). 

For the unbound states with positive energy of the relative motion, E = Jl^q^ /2fj, , scattering phase shifts can 
be found in a similar fashion. For small values of the momentum, qa <C 1 , all phase shifts for m ^ are small 
compared to (s-wave scattering) fP^. The tangent of the s-wave scattering phase shift for q <C ko (i-c E <^Vo) 
is given by : 

t^^So= _ ^ . (12) 



In (cga) 



Jo(Koa) 



Substituting KoaJi{Koa) / Jo{KQa) from Eq. (11) into Eq. (|12D we get 



Note that this expression does not contain parameters of the potential Vq and a explicitly, and it is valid for an 
arbitrary potential well with a shallow m = level [19| . 

Since the integrand in the partition function, Eq. (n^ contains an exponential factor exp{—q'^/q^) , the wavevectors 
q which are larger than the thermal wave vector qt give negligible contribution to the value of the integral. Therefore 
for the short-range interaction or for low temperature, satisfying condition fc^T <C /2^a? , the scattering phase 
shifts need only be considered for q ^ 1/a . Then all the terms in the phase shift sum in Eq. (|l|) can be neglected 
except for the term with m = 0. Finding the derivative of 6q from Eq. (^) we obtain for the two-body interaction 
part of the partition fimction 



cxp(-£;/fcsT) dE 
TT^ + \n^{E/Eb)~E 

The integral in Eq. (|l^) is the Ramanujan integral which can be rewritten as [ pO| , pl] 

dt 



n^ + {\ntY t 



i/(x), (15) 



where 



with X — Eb/ksT . Thus, the partition function acquires a very simple form 

Z^nt = iy{Eb/kBT) . (17) 

A similar result has been obtained recently for contact-interacting particles jw) . 

To consider the small x asymptotic of the function u{x) it is convenient to expand the integral Eq. (|l6|) in 
descending powers of ln(l/a;): 

From Eq. ([l|) one can see that Zint and hence the second virial coefficient B = —X^Zint both vanish when 
Eb/ksT , although one bound state always exists. So the lowest-order density correction to the 2D ideal gas law 
vanishes only slowly as 1/ \n{kBT / Eb) as the binding energy is reduced. Note that, when the potential supports 
several bound states, the contribution of any shallow bound state with m = is cancelled by the scattering phase 
shift integral in the "logarithmic" manner described above. For to ^ the cancellation has a power-law dependence 
in Eij/ksT [ p^ . This implies that higher-order Levinson's theorems responsible for continuity of the partition 
function derivatives p3| are different for m — and m 7^ 0, whereas the zeroth-order Levinson's theorem in two 
dimensions has the same form, Eq. (^), for all to. 

For extremely weak interaction potential, such that (kqa)^ ^ V ln(?i^/2/ia^fcBT), from Eqs. ( |Tl| ) and ( |T^ ) it follows 
that Zint ~ Vofj,a^ /2h'^ , which coincides with the perturbation theory result. In the other limit for large x values 
{Eb/kBT ^ 1), v(x) [Q, therefore the exponential dependence of the partition function on the binding energy 

is recovered. 
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In Fig. |l| we plot the ratio of the total partition function Zint to its bound-state part Zbound = exp{E},/kBT) as 
a function of Ep/kBT. We do this for both the full expression, Eq. (^6|), and the first two terms in the asymptotic 
expansion, Eq. (]l^) . One can see that the asymptotic expression (dashed line in Fig. |l|) is accurate only for very small 
values of Eb/kBT. Over a wide range of E^/kBT both scattering and bound-state terms are important, e.g. when 
Eh = kBT the scattering term produces a 20% correction to When Et/kBT > 3 the bound-state contribution 

dominates completely. 




E/kJ 



FIG. 1. Two-dimensional Boltzmann gas with a short-range weak attraction: the ratio of the two-body partition function 
to its bound-state part, Zi„t /Ztoimd ~ i^{Eb/ ksT)/ exp{Ei,/kBT), is plotted Yersus Et/kBT. Dashed line: the function ^{x) 
is approximated by l/ln(l/2;) + 7/ln'^(l/a:). 

It is instructive to compare the contributions to the second virial coefficient of the direct and the exchange interac- 
tions. For example for a gas of spinless bosons the second virial coefficient is B = — A^(l/4-f 2Zint)- For Zint given 
by Eq. (pT[), the direct interaction term (2Zi„t) is smaller than the exchange term (1/4) for small binding energies. 
Thus, as required for vanishing binding energy, the second virial coefficient is governed by the statistics of particles, 
despite the existence of a single bound state. 
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IV. ELECTRON-HOLE PLASMA 



In this section we study a more realistic model of an interacting gas than the one considered in Sect. [II. We 
consider a mixture of the two types of Fermi particles, positively charged holes (h) and negatively charged electrons 
(e) . This model is important for understanding the thermodynamic and optical properties of the electron-hole plasma 
in semiconductor quantum wells. In the low-temperature, low-density limit most electron-hole pairs are bound into 
excitons. This limit has been studied extensively because of the recent proposals of exciton condensates in this system 
p5| , p6t . The properties of the degenerate 2D electron- hole plasma are also well known . We consider the case when 
the temperature is comparable to the exciton binding energy so the occupation of continuum states (and therefore 
screening) is significant, although the carrier density is low enough to neglect the terms higher than Bn in the 
virial expansion. The temperature and density conditions under investigation are close to those in the regime where 
excitonic gain in wide-gap semiconductors is anticipated [|8| . 

We assume for simplicity equal masses Me = Mh = M and spins s = s' = 1/2 for both species. Then Zhh = Z^e 
and the second virial coefficient for the mixture Eq. acquires the form: 

B = Y Q - 2Ze^ - Ze)j . (19) 

Here A = = {2TTh^ /MkBT)^^'^ and the charge neutrality condition rig = is taken into account. 

The screened Coulomb attraction between electrons and holes is modeled using the Fourier representation of the 
interaction potential: 

27r 

V, = -2^, (20) 

where is the screening wave vector. Hereafter we employ 3D excitonic Rydberg units where length and energy are 
scaled, respectively, by the effective Bohr radius a* and Rydberg Ry* . For electron-electron and hole- hole repulsion 
the same potential with the opposite sign is used. Equation (^) is the well-known Thomas-Fermi expression for the 
Coulomb potential statically screened by a 2D electron gas. Being the long- wavelength static limit of the random 
phase approximation, Eq. ( pO| ) is a simple model for the screened Coulomb potential in two dimensions. Nevertheless, 
this expression reflects the fact that the statically screened potential in two dimensions decreases at large distances 
slower than in the 3D case (as a power law rather than exponentially). Despite numerous more realistic corrections 
|p|,p9| Eq. (^ ) remains the most widely used approximation for the 2D screening |3^,^^. This potential has been 
known for more than three decades [|| but some of its unusual properties were only discovered recently, e.g. the 
existence of a remarkably simple relation between Qs and the number of bound states. Namely, with decreasing 
screening, new bound states appear at the critical integer values of the screening length given by the simple formula 

n ^ (2H+.-i)(2H+^)^ ^^^^2^ _ (21) 

where m is the angular momentum and {v — 1) indicates the number of nonzero nodes of the radial wave function. 



Equation (21 ) shows that several bound states corresponding to the given value of (2|m| -I- v) appear simultaneously. 
This degeneracy is different from the degeneracy for the bound states of the unscreened 2D exciton (or hydrogen atom), 
for which the states with the same value of (|to| -I- j/) are degenerate [^,^. The hidden symmetry that underlines 
this new degeneracy has not been fully understood yet, and Eq. (^ij) still lacks a rigorous analytical derivation. 

We also consider the low-density (nondegenerate) limit, when there is no Pauli blocking and the self-energy correction 
p3| to the Beth-Uhlenbeck formula can be neglected and Eqs. (||j|) and Eq. ^ can be used for Zeh and Zee 
respectively. The shortcomings of this model for the quantitative description of a real system of photoexcited electrons 
and holes in semiconductor quantum wells are self-evident; however, it does provide a tractable model containing all 
the salient features of the system. 

To find the second virial coefficient given by Eq. ( p^ one must calculate the binding energies and scattering phase 
shifts entering the partition functions Zeh and Zee ■ We use for this purpose the 2D modification of the variable- 
phase method [ p4[ known from scattering theory. In this method the scattering phase shift and the function defining 
bound-state energies can be obtained as a large distance limit of the phase function, which satisfies the first-order, 
nonlinear Riccati equation originating from the radial Schrodinger equation. The variable-phase method application 
to scattering and bound states in the screened Coulomb potential ( po[ ) is described in detail in Ref. The method 
is especially effective for calculation of shallow-state binding energies and low-energy scattering phase shifts. 
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Figure |2| shows the results from the calculation of the electron- hole part of the partition function, Ze^,, which 
contains both the bound state sum and the scattering phase shift integral. In this figure Z^h is plotted as a function 
of the inverse screening wave number l/^s- 



40 




FIG. 2. The electron-hole part of the partition function, Z^h versus the screening length 1/qa for two values of fcsT/Ry. 
Solid lines show the bound state contributions Ztound only. Dashed lines: Ztound + Zscatt- 

For the electron-hole plasma is a function of carrier density |3C|], and in the purely 2D case Qs oc (jig -I- 12^) 
for low densities and is independent of density, Qs — > 8/a* for Me = Mh, in the degenerate limit. Here we treat qs 
as a parameter characterizing the strength of the screened interaction potential Eq. To emphasize the role of 

scattering we show on the same plot the bound-state sum, Zbound = X^m j/ 6xp(— i^^.i/Z^B?^) , which exhibits jumps 
whenever new bound states appear [i.e. when Qs satisfies Eq. (pi])]. These jumps become higher with increasing 
screening length I/qs since several bound states appear simultaneously. As can be seen the additional scattering 
phase shift contributions completely remove these jumps. The partition function is plotted for two values of the ratio 
of ksT to the excitonic Rydberg, ksT — 1 Ry* and ksT = h Ry*, which roughly correspond to ZnSe (or GaN) 
and GaAs at room temperature. One can see that for high temperature (or low binding energy) the bound-state 
contributions to Z^h are suppressed by the scattering phase shift integral more strongly than in the fc^T = 1 Ry* 
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case. 

In Fig. ^ the second virial coefRcient B (scaled by A^) is plotted versus the screening wave number qs for two 
different values of fc^T. 




FIG. 3. The second virial coefficient of the dilute electron-hole plasma B/X^ as a function of the screening wavenumber 
Qs ■ Solid line: ksT = 1 Ry*; dashed line ksT = 5 Ry*; dot-dashed line: non-interacting dilute electron-hole plasma, 
B/X^ = 1/16. 

Equation ( p^ is used for the calculation of B, and the repulsion (Zee) term partially compensates the Z^h term. 
This compensation is especially significant in the high-temperature case ksT = 5Ry*, in which the 2D electron- hole 
plasma behaves much like an ideal gas over a wide range of screening wave vectors. For ksT = 1 Ry* the electron-hole 
attraction term dominates and the plasma is strongly correlated for all values of qs- In this case a small statistical 
repulsion {B/X^ = 1/16, horizontal line in Fig. ^), which is due to the fermionic nature of electrons and holes, can 
be completely neglected. Thus, at room temperature the electron-hole plasma in GaAs-based quantum wells can be 
treated as an ideal gas, whereas in wide-gap semiconductors (e.g., ZnSe or GaN) due to the high value of Ry* the 2D 
electron-hole plasma is strongly correlated and excitonic effects are important for its thermodynamic properties. 
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V. CONCLUSION 



In this paper we show that a proper account of scattering ehminates discontinuities in thermodynamic properties of 
the nonideal 2D gas whenever extra bound states appear with small increase of the strength of an attractive potential. 
This treatment provides a heuristic proof of Levinson's theorem in two dimensions. 

We trace the way in which the bound-state contribution to the partition function of the 2D gas, with a weak 
short-range attraction between its particles, vanishes when the binding energy decreases. A weak \/ IniksT / Eb) 
binding energy dependence of the second virial coefficient of such a gas is found for Eh/kBT — > 0. 

A consistent treatment of bound and scattering states in a screened Coulomb potential allows us to calculate 
the quantum-mechanical second virial coefficient of the dilute 2D electron-hole plasma and to establish the difference 
between the nearly ideal electron-hole gas in GaAs and the strongly correlated exciton/free-carrier plasma in wide-gap 
semiconductors . 

The 2D electron-hole plasma was considered in the low-density nondegenerate limit only. Transition to the strongly 
degenerate Fermi limit and related questions of Pauli blocking and self-energy corrections to the Beth-Uhlenbeck 
formula in 2D remain the subject of further research. 
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